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A  simple 

transfer  rate  from! 
tion  6(  fluid  propej 
behind  a  reflected! 
in  thermodynamic! 


^  SUMMARY 

P  '  ' 

111  method  is  developed  for  calculation  of  the  heat 

i 

Uitrlty  boundary  layer  equations  with  arbitrary  varia- 

fe'  '  ' 

pThe  method  is  applied  to  the  end  wall  of  a  shock  tube 
||ln  a  perfect  gas,  and  to  the  stagnation  point  for  a  gas 
Ibrium  with  Lewis  number  unity.  Explicit  formulas 


and  results  are  oblwSpA. when  the  fluid  properties  are  powers  of  the  enthalpy 


or  temperature,  e| 
For  the  end 

I 

fer  rate  is  derived| 
3%  for  all  caset^ 
'For  the 

equations  are  obt^ 
the  momentum  Ian 


is  thicker,  Numei 


correlated  with  e: 
the  pressure  graj 
only  one  parameti 
this  correlation  ^ 


A  Prandtl  number  is  constant. 

P' 

1  geometry,  an  analytical  expression  for  the  heat  trans- 
sb  can  be  correlated  with  exact  calculations  within 


Mon  point  geometry,  two  pair  of  transcendental  algebraic 
K  one  pair  of  which  gives  the  heat  transfer  rate  when 
tf  thinner  than  the  energy  layer,  the  other  pair  when  it 


solutions  of  these  equation 


can  be 


(elutions  to  within  ^  4%.  The  integral  method  indicates 
arameter  and  Prandtl  number  can  be  combined  into 
d  examination  of  results  of  exact  calculations  shows 
lite  accurate.-—^ 
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NOTATION 


constant  of  integration,  Eq.  (3.20) 
specific  heat  per  unit  mass 

frozen  specific  heat  per  unit  mass  of  gas  mixture 
dimensionless  stream  function,  Eqt*  (3*  8) 
dimensionless  total  enthalpy,  Eq.  (3.8) 
total  enthalpy 

integral  expression,  Eq.  (3.25) 

Eq.  (3.40) 

Eq.  (3.  42) 

index  of  dimensionality,  =  1  for  eocisymmetric  flow 

=  0  for  two-dimensional  flow 


K  pk 

k  thermal  conductivity  (frozen,  in  the  case  of  an  equilibrium  gas) 

N  distance  from  edge  of  energy  boundary  layer,  Eq.  (3.31) 

p  pressure 

-q^  heat  transfer  rate  to  surface 
T  absolute  temperature 

t  time 


u  velocity  component  parallel  to  surface 

V  velocity  component  normal  to  surface 

X  coordinate  parallel  to  surface 

y  coordinate  normal  to  surface 

q  similarity  variable,  Eqs.  (2.7)  and  (3.  II) 


! 

•  i 
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.  SECTION  1 
INTRODUCTION  ^  ^ 


Many  pf  the  heat  transfer  boundary  layer  problems  of  interest  in  con- 
nection  vyith  high-speed  missile  and  space  vehicle  flight,  as  well  as  related 
laboratory  cohiiguratibns,  involve  similarity  solutions  of  the  boundary  layer 
equations^  Exarnples  are  the  stagnation  point,  the  end  wall  of  a  shock  tube, 
and  the  coristarit-presSure  flows  along  cones,  wedges,  and  side  walls  of  shock 
tubes.  Similarity  solutions  are  also  important  as  the  building  bldcks  of  the 
"Ibcal-similarity"  method  of  approximate  boundary  layer  computation,  which 
has  Come  into  wide  use  in  recent  years:  : 

The  usefulness  of  the  similarity  equations  is ’as sociated  with  the  fact 
that  they  represent  a  reduction  of  partial  differential  equations  in  two  varia-  . 
bles  to  ordinary  differential  equations,  thus  facilitating  immensely  the  numer¬ 
ical  work  needed  to  obtain  solutions.  Another  method  of  performing  this 
reduction  to  ordinary  differential  equations  has  long  been  used  in  boundary 
layer  theory  --  the  Karman-Pohlhausen  method.  In  this  technique,  the  par¬ 
tial  differential  equations  are  integrated  once  in  a  direction  normal  to  the 
surface,  leaving  ordinary  differeritiaTequatioris  involving  integrals  of  the  de-.  , 
pendent  variables  such  as  velocity,  temperature,  etc.  Profiles  for  these 
quantities  are  assumed  and  inserted  in  the  integrals,  enabling  evaluation  of 
the  coefficients  of  the  ordinary  differential  equations.  . 

The  idea  of  this  paper  is  to  combine  the  two  techniques,  treating  the 
ordinary  differential  equations  of  similarity  theory  by  the  Pdhlhausen  integral 
method,  and  so  arriving  at  an  expression  for  the  heat  transfer  rate  to  the 
surface  involving  only  simple  quadratures.  This  application  of  integral  meth- 
qds  to  ordinary  differential  equations  has  been  used  before  in  some  special 
cases.  The  incompres sible  momentum  equation  for  the  flat  plate  and  the  two- 
dimensional  stagnation  point  are  discussed  in  Schiichting' s  fextbook.  V  Bush^ 
uses  it  on  the  similar  compressible  boundary  layer  equations  for  the  case 
where  the  product  of  density  and  viscosity  is  constant  and  the  Pfandtl  number 
is  unity.  It  is  our  object  here  to  develop  the  method  for  arbitrary  fluid  prop¬ 
erty  variations  (with  examples  worked  Out  for  power  law  variations)  with  an 
eye  toward  later  application  to  heat  transfer  problems  involving  dissociated 
and  ionized  gases. 

The  basis  of  the  method  is  the  representation  of  the  temperature  (or 
enthalpy)  profile  as  a  function  of  the  similarity  variable  by  an  expression 
which  is  simple  enough  to  be  easily  integrable,  yet  which  leads  to  values  for 
the  heat  transfer  which  are  sufficiently  accurate.  .  The  profile  used  here  was 
first  suggested  byTepson,  ^  who  obtained  a  remarkably  accurate  formula  for 
heat  transfer  to  a  flat  plate  at  Prandtl  number  zero,  with  the  thermal  con¬ 
ductivity  an  arbitrary  power  of  the  temperature.  It  was  this  v/ork  which 


stimulated  the  developments  described  in  the  present  report. 

Two  cases  will  be  discussed  to  illustrate  the  method  of  solutioh 
the  similarity  equations.  The  first  is  the  heat  transfer  frdna  a  perfect  gas 
to  the  end  wall  of  a  shock  tube;  the  second  is  the  heat  transfer  from  a  gas 
in  therrriodyhamic  equilibrium  with  Lewis  nurriber  unity  to  a  stagnatioh'point. 
When  power  law  fluid  properties  are  Used,  the  first  case  leads  to  a  simple 
algebraic  forrhuia  for  heat  transfer  irate;  and  the  second  leads  tu  a  pair  of  . 
franscendental  algebraic  equations  that  can  be  Solved  rnuch  more  quickly 
than  the  similarity  differential  equatiOhs,  which  require  the  ihtegratipn  of 
a  two-point  boundary,  value  problern.  ^  . 

There  is  no  treatment  in  this  report  of  cases  where  the  energy  equa¬ 
tion  is  coupled  to  continuity  equations  for  various  species  in  the  gas,  as  would 
be  the  case  in  a  dissociating  or  ionizing  gas  which  did  hot  fulfill  the  conditions 
of  thermodynamic  equilibrium  and  Lewis  number  unity.  However,  extension 
of  the  method  to  this  situation  is  under  investigation  and  will  be  the  subject 
of  future  reports.  ■ 

The  author  r.vould  like  to  acknowledge  the  frequent  advice  and  encour¬ 
agement  of  Prof.  James  A.  Fay  during  the  course  of  the  work  reported  here. 
The  machine  computations  were  carried  out  by  Mr.  Richard  Goodman  and 
Miss  Patricia  Jackson. 


SECTION  2 

END  WALL  HEAT  TRA.NSFER 


Between  the  end  wall  of  a  shock  tube  and  the  reflected  shock  which 
recedes  from  it  there  is  a  sample  of  hot  gas  which,  neglectihg  boundary  layer 
displacernent  effect,  is  stationary.  Being  in  contact  with  the  cold  wall,  the 
gas  is  cooled  by  a  tempe ratur e  boundary  layer  g rowing  from  the  wall.  This 
is  a  Rayleigh-type  problem,  since  it  can  be  represented  by  a  cold  wall  sud¬ 
denly  put  in  contact  with  a  hot,  quiescent  gas  sample.  If  the  gas  \yere  treat¬ 
ed  as  incompressible,  we  would  expect  the  solution  to  depend  on  the  usual 
similarity  variable  of  Rayleigh  problems,  y/./T,  where  y  is  distance  normal 
to  the  wall-  Because  of  compressibility,  we  expect  that  this  variable  is 
modified  by  the  Howarth  transformation  to  /  p  dy/ -JT. 

Let  us  first  consider  a  perfect  gas  with  no  dissoCiafive  or  ionizing 
effects.  There  is  no  flow  parallel  to  the  wall(in  the  x  direction)  nor  are  there 
any  derivatives  in  this  direction.  The  continuity,  normal  momentum  and 
energy  equations  are 


dt 


^  o 


{2-  la) 


(Z.lb) 


(2-  Ic) 


The  continuity  equation  indicates  that  density  changes  lead  to  velocities  nor 7 
mal  to  the  wall.  When  the  usual  boundary  layer  order  of  magnitude  arguments 
are  applied  to  these  equations  (2.  lb)  tells  us  that  p  ”  Pg  ^  )  where  pg  is 

the  pressure  external  to  the  boundary  layer,  so  that  the  pressure  terms  in 
(2.1c)  can  be  ignored.  The  dissipation  term  in  (2.  Ic)  is  of  order  p  =  0(62) 
also,  and  thus  also  negligible.  The  boundary  layer  energy  equation  then 
becomes--'-' "■ 


This  combined  with  the  continuity  equation  (2.  la)  and  an  equation  of  state 
serve  to  define  t^eiemperature  and  flow  field.  The  boundary  cdnditm^  3-re: 


V.0,  1-  =  Vc 


oris  tan  t 


{2.3a) 


T=X  =  constant 


{2.3b) 


The  velocity  V  can  be  eliminated  by  use  of  the  continuity  equation 
(2.  la).  ^This  enables  us  to  write  the  convective  derivative  operator  as 


We  now  introduce  the  similarity  variable  by 


{2.  4) 


{2.5) 


and  assume  that  T  arid  p  depend  on  q  alone.  If  q  is  then  inserted  into  {2.  4) 
and  {2.  2),  the  result  is 


A  f  li  ill  — i-  — 

‘‘’i  <*>1  '  V  X|  J-t  X,  f.fe. 


=  0 


{2.6) 


where  the  subscript  o  refers  to  quantities  evaluated  at  same  reference  condi¬ 
tion -to  be  chosen  later.  The  quantity  Xj  is  chosen  to  make  the  expression  in 
parenthesis  in  the  second  term  equal  to  Unity,  yielding 


-Vzri 


{2.7) 


If  we  further  non-dimensionalize  T  with  the  external  temperature  T^,  the 
final  equation  we  must  solve  is 


i.  (b  4®!  +  Al  t]  P-O:  T  ' 


{2.8) 


with  boundary  conditions 


,  ,  e(cD)=i  :  •  ‘ 

The  heat  transfer  rate  to  the  wall  can  be  written  \vith  the  helo  of 


We  see  that  a  heat  flux  potential  p ^ ,  may  be  defined  by 

.  /  &  ■ 

d<i>,  _  K  ae  i-  i  ,  f  Ji 


such  that  the  heat  transfer  rate  to  the  wall  is 


and  the  derivate  of  pj  replaces  the  derivative  of  T  for  calculating  q^. 

The  introduction  of  a  heat  flux  potential  is  a  frequent  device  in  heaV 
transfer  problems,  and  has  been  applied  to  gases  in  the  end  wall  geometry 
in  Ref.  4.  However,  in  this  reference  convection  was  neglected,  which  is 
equivalent  to  ignoring  compressibility,  so  the  density  factor  did  not  appear 
and  p  became  the  classical  potential  /  kdT.  It  would  appear  that  the  omis¬ 
sion  of  convection  casts  doubt  on  the  validity  of  the  work  presented  in  Ref.  4. 


We  now  propose  to  solve  Eq.  {2.  8)  by  an  integral  method.  To  this 
end  we  integrate  on  q  from  q  =  0  to  q  =  «>,  remembering  that  de/dq  —  0  as 
q—*  =0.  The  result  is 


(2.  13) 


This  form  shows  that  all  we  need  to  find  is  an  assumption  for  Cp  and  q 
as  functions  of  0.  Notice  that  0  is  now  the  independent  variable  in  the  in¬ 
tegral,  and  the  range  of  integration  is  finite.  This  fact  is  of  crucial  inipof- 
tance  in  the  approximation  scheme.  The  relation  between  c^  and  O  (or  T) 
is  given  by  thermodynamics.  The  chief  question  is  \vhat  to  use  for  q  (0}; 
For  that  ve  return  to  the  differential  equation  (2:  8)  andhotice  that  near  the 
■wall''whege'-q 'if- '^.imall  '■ 


^  A  —  n 


(2.  14) 


W e  nOw  chobse  the  q  profile  by  integrating  this  relation  to  give 


{.Ill 


(2.15) 


This  is  not  a  linear  relation  between  q  and  O,  but  a  linear  relation  behyeen 
q  and  the  heat  flux  potential,  which  we  will  see  gives  a  power  relation  between 
q  and  0  for  the  case  where  k  is  a  power  of  0.  It  was  for  that  Case  that  the 
profile  (2,  15)  was  introduced  by  Jepson.  ^ 

By  inserting  (2.  15)  in  (2. 13),;  using  the  definition  of  in  (2.  11),  we 
find,  the  final  form  ' 


Ko  / 


(2.  I6a) 


This  is  a  simple  quadrature  formula  for  the  heat  transfer  rate.  In  the  case 
where  Cp  may  be  taken  constant  integration  by  parts  leads  to  the  simpler 
expression 


(2. 16b) 


To  summarize  the  procedure,  we  integrate  as  in  the  ordinary 
Pohlhausen  method,  but  use  the  dependent  variable  O  as  the  variable  of 
ihtegfation.  We  must  then  insert  a  profile  q{0).  We  derive  the  profile  by 
evaluating  the  differential  equation  near  q  =  0,  and  integrating  that  expres¬ 
sion.  We  note  that  the  profile  so  derived  does  not  satisfy  the  condition  ; 


-6- 


as  which  was  . juiiTiined  in  the  derivation  of  the  inte:gral  condition 

(2.  13).  However,  the  protileris  only  used  in  the  intogrhl  bn  0,  and  the  error 
that  arises  in  only  caused  by  the  difference  between  the  areh  under  the  exact 
^(0)  curve  and  the  (pj  -  vs  0  curve  between  tlie  finite  limits 

Q  =  0^^  and  0  =  I,  Assessment  ot  the  actual  error  is  Obfained  by  Gomparisbn 
of  Specific  cases  with  exact  solutions. 

Application  to  Power  Law  Thern'.al  Conductivity 


To  investigate  the  accuracy  of  the  approximate  formula  (2.  I6b)  we 
turn  to  several  specific  cases.  For  a  perfect  gas,  Cp  is  constant  and  the 
equation  of  state  in  the  present  constant  pressure  situation  indicates  that 
p  is  inversely  proportional  to  T,  Let  us  take  k  proportional  to  T'T ,  a  fre¬ 
quent  representation  of  thermal  conductivity.  For  convenience  \ve  also  , 
\vill  choose  the  reference  condition  o  to  be  the  external  condition  e.  Then 
Eq.  (2,  8)  becomes 


V 

r^-e 

=  0 

and  from  (2,  11)  and  (2.  l6b) 


(2.  18) 


^V= 


I”  _  l~ 

t*j  OM-t 


U.  19) 


This  is  the  result  obtained  by  Jepson.  ^  .  Notice  that  pj— “0"  ,  so  the  assumed 
relation  between  ri  and  0  is  a  linear  one  for  q  in  the  terms  of  0'^  ,  as  rrien- 
tioned  above. 

The  formula  (2.  19)  yields  a  finite  heat  transfer  rate  as  ©vv  T  0  even 
for  <  1,  which  would  not  be  true  if  a  linear  relation  between  q  and  ©  had 
been  used  for  the  profile.  Another  important  limiting  case  is  io  =  1,  which 
corresponds  to  pk  constant,  for  which  we  have 


(2.  20) 


This  same  expression  is  Obtained  in  the  limit  i-  In  both  these  cases 

ah  exact  sblutioh  of  Eq.  (2.  17)  is  possible,  and  yields 


A  comparison  of  formula  (2.  19)  With  exact  solUtiohs  of  Eq.  (2,  17)  is  given  in 
Table  1/  The  exact  calculations  are  from  Jepson^  and  from  unpublished  re¬ 
sults  obtained  in  the  course  of  an  investigation  by  Adams, 


TABEE  I 


■ 

e 

w 

Appro.ximate  {dpi/dq)^^ 
from  (2.  19) 

Exact 

(dpi/dn)w 

;Ref;^ 

1/2 

■  0 

,  _  1.154 

1.297 

■  ,  l.:i2'-f--- 

1/2 

.  0375 

;  .  975  : 

l.OSO 

11 

1/2 

•1 

.850 

.950 

1.  12 

1/2 

.395 

.442 

1.  12 

3 

1/2 

.  75 

■  :  .185'  ; 

.  208 

1.  12 

,  ,3;:' 

5/2 

0 

^.338'' 

V.  :  .389,;,:, 

'  ■■■  :■  1-  15  ' 

5/2 

.0375 

.338 

.  390  . 

1.  16 

5/2 

.1 

•  336 

.  386 

1.15 

3 

5/2 

.85 

.319 

.371 

1.  16 

■  ,'3_.: 

5/2 

.5 

.  . 262' - 

.298 

1.  14 

■',3' 

5/2 

175 

.  154 

.176 

1.  13 

3 

1 

©w 

(1  -  g,^)/vT 

(1  - 

=  E  13 

.—  I-  ■ 

(1  -  0,^)  /  ^ 

(1  - 

-v/4/-  ,=  -.,1,13  ■ 

We  see  that  in  all  eases  (2,  19)  yields  values  smaller  than  the  exact  values, 
but  the  ratio  of  the  exact  to  the  approximate  solution  is  remarkably  constant 
over  a  wide  range  of  w  and  0^. 


Another  case  for  which  exact  calculations  are  available  is  one  in 
which  the  k  vs  T  curve  is  represented  as  consisting  of  two  different  pov/er 
laws  with  a  break  at  some  temperature  T^,.  Adams^  used  this  form  in  his 
study  of  electronic  heat  transfer.  The  relation  for  the  k  curve  is 

lor 

=  ©  for 


This  leads  to 

r  ^  v-x  x-i .  f 

A  C^-e)  et  e  a©  v  1 

Vci'lV 

t 

0 

% 


v-\ 


die 


X*l 


V  .  .  ■ 

1 


(2.  23) 


+A-  [vCi^+l)  ^ 


S^\ 

— 1 
X4 1  j 


This  of  course,  has  the  same  limit  for  0.,,  1  as  the  single  exponent  case, 

■  Eq.;:;(2;21). v:.' 

Comparison  of  Eq.  (2.  23)  with  unpublished  exact  calculations  of 
Adams  are  shown  in  Table  II.  . 

■  TABLE  II"' 


■ 

■ 

0 

w 

0 

c 

Approximate  (dpj/dq)^ 
from  (2.  22) 

Exact 

(dpi/ dTi)^v 

Exact 

Approx. 

1/2 

5/  2 

.222 

.  889 

.Vv  .602,  ■ 

.  67  5 

■''i.w'W: 

1/2 

5/2 

.167 

.667 

509  ■  ■  ■ 

.  580  . 

M.. 14 

1/2 

5/2 

.125 

.500 

.459 

.  505 

1. 15". 

X 

V 

—1 

—  1 

(1  -  e„)/VI 

V2A  (1  -  0  ) 

w 

•/lA  =  1. 13 

Again  we  see  the  approximate  solution  (2.  23)  gives  smaller  values  than  the 
exact  one,  but  a.gain  the  ratios  are  very  constant. 

On  the  basis  of  the  evidence  presented  in  Tables  I  and  II.  it  appears 
that  multiplication  by  the  factor  A/ -l/“,  which  makes  the  appro.ximate  solution 
agree  with  the  exact  one  in  the  Rayleigh  case  0^ -*■  1,  gives  an  answer  correct 
to  within  +  3%  for  all  the  numerical  cases  calculated.  Therefore  the  corrected 
formula  (T.  12)  for  the  heat  transfer  rate  becomes  (with  external  reference 
-conditions). 


2.t 


X, 


I 


o,., 


(2.24) 


This  provides  a  simple  and  accurate  formula  for  Calculating  the  heat  transfer 
rate  to  the  end  wall  of  a  shock  tube  when  p  and  k  are  given  function  of  T,  and 
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SECTION  3 

STAGNATION  POINT  HEAT  TRANSFER 


Heat  transfer  to  the  stagnation  point  of  a  bliint  body  has  been  studied 
intensively  in  recent  years.  It  also  represents  a  case  which  is  exactly  re¬ 
ducible  to  similarity  form.  Here  we  shall  restrict  Ourselves  to  the  case  of 
a  dissociating  gas  with  Lewis  number  unity  and  flow  in  thermoSynarhic 
equilibrium,  to  illustrate  the  application  Of  the  method  withbut  undue  compli 
cations.  Under  these  circumstances,  we  have  to  solve  the  following  contin¬ 
uity,  mbm.entum,  and  energy  equations^  for  the  velocities  u  and  v  and  the 
total  enthalpy'H:'  ; 


=  o 


(3.  1) 


2X 


(3.2) 


■jh  ,'3H 


(3.  3) 


k  and  are  the  frozen  thermal  conductivity  and  specific  heat  of  the  gas. 
The  dissipation  and  pressure  gradient  terms  do  not  appear  in  the  energy 
equation  (3.  3)  because  they  are  negligible  at  a  stagnation  point.  The  factor 
rJ  with  j  =  0  for  two-dimensional  and  j  =  1  for  axi- symmetric  flow  permits 
consideration  of  both  cases  at  once,  with  r  the  cylindrical  radius  of  the  axi- 
symmetric  body.  Near  the  stagnation  point  r  «  x.  These  equatioris  define 
the  problem  when  combined  with  the  boundary  conditions 


r-O 


|=0D 


H=Hw  - 


a=u. 


constant 


(3.4) 


(3.5) 


where  near  the  Stagnation  point  we  have  put  r  =  x,  ug  =  x{due/dx)st- 
this  choice  of  \2  >  and  putting  K  s  pk/cp,  the  similarity  equations  (3.  9)  and 
(X  10)  becorfie 


I 


(3.  12) 


(3.  13) 


while  the  boundary  conditions  (3.  4)  and  (3.  5)  transform  to 

>^=o: 


r  00 


(3.14) 
(3.  15) 


Equations  (3.  12)  and  (3.  13)  differ  slightly  from  the  usual  stagnation- 
point  boundary  layer  equations,  ^  in  that  the  Prandtl  number  or  =  Cp 
appears  in  the  momentum  equation  instead  of  the  energy  equation.  This  is  a 
consequence  of  choosing  the  quantity  K  =  pk/Cp  as  the  basic  combination  of 
fluid  properties,  rather  than  the  more  usual  choice  of  pu.  Since  the  heat 
transfer  rather  than  the  shear  stress  is  of  primary  interest,  it  seems  more 
suitable  to  relegate  the  Prandtl  number  to  the  momentum  equation,  which  is 
only  weakly  coupled  to  the  energy  equation.  The  use  of  pp  is  a  holdover  from 
the  point  of  view  that  the  momentum  equation  was  basic  and  the  temperature 
profile  could  be  found  after  the  velocity  profile  was  known.  For  the  case  of 
constant  Prandtl  number,  Oq  .  (3.  12)  and  (3.  13)  can  be  put  into  their  more 
familiar  form  by  replacing  q  by  .y!^q  and  f  by.y^f. 

The  heat  transfer  rate  is  expressed  by 


(3.16) 


In  this  case  we  may  again  define  a  heat  flux  potential  A  by  inserting  ^2 
from  (3.  il)  into  (3.  16): 


(l+i)  fo 

(3.17) 
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The  problem  before  MS,  of  determining  {d(J2/dT|)^  from  Eqs.  (3.  12) 
and  (3, 13V  plus  some  expressions  for  the  fluid  properties  in  equilibrium  as  a 
function  of  g,  is  more  difficult  than  in  the  end  wall  gebmetry  because  of  the 
presence  in  the  energy  equation  of  the  stream  function,  f,  determined  by  the 
tnomehtum  equation.  A  proper  solution  requires  simultanebuS  conSideratibn 
of  both' equations.  Howev’Cr,  it  is  well  known  that  the  influence  bf  the  tnbmen 
turn  equation  bn  the  energy  equation  is  small;  arid  so  we  may  hope  to  include 
this  influence  by  an  approximate  treatmeiit  of  Eq.  (3.  12).  This  will  be  the 
-  appr bach"  we ‘Shall 'take'. . 


First  consider  the  energy  equation  (3,  13)  and  integrate  from  the  wall 
to  the  edge  of  the  thermal  boundary  layer,  where  dg/dq  =  O.  Using  the  def- 
initiorr  (3. 17)  of  and  integrating  by  parts,  we  find 


an  equation  analogous  to  {2.  13).  We  must  now  have  a  guess  for  the  profiles 
df{g)/dq  and  n  (g)  just  as  in  the  end  wall  case  we  needed  one  for  q  (0  ).  The 
momentum  and  energy  equations  will  supply  these  profiles. 


The  f  profile  is  determined  much  in  the  iyay  we  determined  the  0 
profile  in  the  end  wall  case,  by  writing  (3.  12)  near  the  wall  as 


i.  (  S  A  ^  O 


(3.  19) 


We  now  integrate  to  find 


d‘+  - 

J_  Xo 

r-PMan  +  P 

(3.20a) 

df  ■ 

L  ^ 

f  1 

1-i 

(3.20b) 

1  y 

L  Jo  l 

J  cr- 

where  A  is  a  constant  of  integration,  related  to  the  shear  stress,  which  will 
be  evaluated  later. 


wall: 


The  t)  {g)  profile  is  obtained  the  same  way  from  Eq.  (3.  13)  near  the 


i-  =0 


Xo 


^  _  f 

mi 


.fik 


(3.21) 


This  leads  to 


(3.22) 


^^ener% 

profile  is  g  =  1  Jnstead^;hl  ;  layer^  the  g  . 

will  have  to  consider  separatelv  the^c'^  °  \  Thus,  we 

than  the  momentum  lave^  and  tL  ®  the  energy  layer  is  thicker 

These  relative  thicknesses  d^n  f  t  momentum  layer  is  thicker, 

and  will  prov^L  one  oYthe  n  r  ’  °  "^"^her  a. 

ties  dep^cTom^  Zy  problem.  Since  the  fluid  proper- 

of  the  energy  layer  is ’spec if ild  a^  T""  variable.  Then  the  edge 

momentum  Uyor  is  the  unknown,  detfned  ai  g  e",  “fi,'”" 

energy  layer  If  it  \<  thir-\r^^  *^u  °  thinner  than  the 

that  fh’;  drata;cl\‘'.we‘.'’r,irete'’3  o^'rrio^ 

venient  variable.  '  °  layers,  nf  -  T|g,  is  a  con- 

eral  meflmd^  solSon  '^e^a^  us  consider  the  gen- 

1-imary  interest  is  the'heat 

depending  onihe^h^^^^^trnTt'o^t;;ick:^;,rnTh^^  eL' ‘'y'LK 
Conditions  are  therefore  needed  for  the  solution.  We  have^ail^; 

*u  1  .  '  .  '  '  '■u  ine  eoge  oi  the  momentum  laver  This:  can  n..  •  j 

through  and  vields  a  raf-Vin,.  inne.t.u  .•  •  ,  yvr.  inis  can  be  carried 

B  q  yields  a  rather  lengthy  equation  involving  A  and  (d^2/dT,) 
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algebraically  and  transcendentally  in  the  lirnits  of  integration.  However, 
a  simpler  approach  to  the  determination  of  A  was  sought,  in  view  of:^he"fact 
that  the  only  appearance  of  f  is  under  the  integral  sign  in  Eq.  (3/ 18)  as 
df/drj  j  so  Only  a  weighted  average  value  of  df/di^  is  computed.  We  have  al¬ 
ready  imposed  oh  df/dri  in  (3.  20)  the  condition  that  it  vahish  at  the  wallj  and 
satisfied  the  itiomehtum  differential  equation  hear  the  wall,  and  we  will  use 
the  condition  that  it  is  unity  at  the  edge  of  the  rhomenturh  layer.  In  lieu  of 
making  it  satisfy  the  rnomentum  equation  On  the  average  throisghout  the  whole 
momentum  layer,  we  instead  impose  a  Second  boundary  cpnditiOn  at  the  edge 
of  the  layer,-' making*,"  ■ 


-  0  h  >1 


(3.23) 


This  serves  as  the  third  relation,  which  completes  the  determination  of  the 
problem,  ;$ince  this  condition  is  largely  instrumental  in  determining  A, 
which  in  turn  is  proportional  to  the  shear  stress,  one  will  not  expect  to  de¬ 
termine  that  quantity  accurately.  However,  since  df/dq  now  has  the  correct 
values  at  T]  =  0  and  q  =  qf,  and  the  correct  derivative  at  t)£,  and  satisfies 
the  correct  equation  hear  ri  =  0,  we  may  hope  that  the  effects  of  the  momen¬ 
tum  layer  on  the  heat  transfer  are  sufficiently  well  represented.  Eventual 
comparison  with  exact  sqlutipns  will  provide  us  with  the  ultimate  check  on 
the  usefulness  of  this  simple  approach.  . 

With  the  use  of  the  third  condition  (3.  23)  established,  we  may  now  go 
on  to  consider  separately  the  cases  of  morhentum  layer  thinner  and  thicker 
than  the  energy  laye r.  "  ■  ■ 

Momentum  Layer  Thinner  Than  Energy  Layer  (Small  Prandtl  Number) 


The  momentum  layer  terminates  at  g  =  gf  <  1.  From  gf  to  g  =  1, 
the  velocity  is  free  stream;  i.  e.  df/dri  =  1.  From  g^  to  gf  Eqs.  (3.20)  are 
applicable.  The  constant  A  is  determined,  as  discussed  above,  by  letting 
d^f/dq^-  O  at  q  =  qf.  The  use  of  this  in  Eq.  (3.  20a)  gives 


(3.24) 


We  now  define  the  integral 


(3.25) 


r This  idea  grew  out  of  a  convers 
Dr.  Richard  Levy. 


tion  with  the  author's  colleague. 
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where  the  second  form  comes  from  nsing  (3.  22).  Then  (3.  2 Ob)  becomes 


where  l£  ^Sf)  ^”d  (3.  22)  has  been  iised  again  to  replace  di^ . 


(3.26) 


The  boundary  condition 


i  1  at  Ti£  (or  g^)  then  leads  to  the  condition 


(3.27) 


The  basic  expression  (3.  18)  for {di?2/dTi ),.v  can  now  be  written  in  final 
form  by  breaking  the  range  of  integration  at  g£  and  usin|  (3.  22)  for  drj/dg 
and  (3.  26)  for  di/dr]  in  the  lower  part  and  (3.  22)  and  df/dri  =  1  in  the  tipper 
part.  The  result  is 


a  41  ■ 

.  r  a  -I 

1+.J  Xo  o-J  / 

If  the  last  term  is  integrated  by  parts  and  (3.  27)  is  used,  this  may  be  written 


Mi/wKdiL  Jc.,  ^  ^ 


(3.  28b) 


Equations  (3.27)  and  (3,28)  are  two  relations  for  the  two  unknowns 
(d02/dr|  add  g£.  Once  K  and  <j  are  given  as  functions  of  g,  they  can  be 
solved  and  the  heat  transfer  rate  determined  from  (3.  17). 

As  the  Prandtl  number  CT  —  0,  the  momentum  layer  thickness  vanishes 
so  gf  -»  g.j^.  The  third  term  of  (3,28b)  disappears,  and  the  equation  reduces 
to  Eq.  (2.  lob)  for  the  end  wall  geometry.  This  is  the  physicaliy  correct 


(^1  -/MC 


r  1 


V  X, 

i.' 


:] 


Applicatiori  of  the  boundary  condition  finally  yields. 


. 


^  r  ^ 


- 


3 


(3.35) 


The  basic  integral  for  (d6>/ dr|)^y,  (3.  18)  is  now  written  iis^tng  (3.  22) 
for  dT|/dg  and  (3.  33)  for  df/dq .  The  result  is 


(3.36) 


The  two  relations  (3.  35)  and  (3.36)  determine  the  unknowns"{36>/dT| )^ 
and  X£,  the  latter  replacing  gf  in  the  present  case.  Again  specification  of 
K  and  a  as  functions  of  g  enable  us  to  solve  for  the  unknowns  and  thus;  the 
heat  transfer  by  , (3.  1 7).  . 

.Applicatioh  to  Power  Law  Fluid  Properties 

In  order  to  test  the  present  appro.x:imate  method,  calculations  have 
been  made  for  constant  Prandtl  number  CTq,  and  the  fluid  properties  have  been 
taken  as  power  functions  of  the  enthalpy,  with  the  reference  condition  as  the 
edge  of  the  energy  boundary  layer,  i.  e.  , 


o>-l 


5  =  3 


(3.  37) 


This  simple  representation  permits  the  integrals  to  be  calculated  analytically. 
The  boundary  layer  equations  (3.  12)  -  (3.  14)  for  this  case  become 


(3.  38) 
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For  the  first  case,  where  the  rnomentum  layer  is  thinner,  Eqs.  (3.  27) 
and  (3.  28a)  become 


I'ii'l  -  1 

_  q 

-X-to+2.  -X  ooi+.£ 

( 


r  i-  jT  _  I- 


\  3  /w  1  to 


«**+! 

34^ 


to-V( 


! 


-1 


no; 


T"  »*>>  i 

a>4.( 


2.10+2 


(3.39a) 


-l‘  V- 

\-^to+jL  -l-«04-2.  J 

ta+^  ■ 

(irh  -  vC') 


(3.  39b) 


where 


h 


Xn^  I  , 


(3.  40) 


To  verify  the  approach  to  the  end  wall  result  as  aQ  —  O*  (Xe.  ,  g£-^  gw)» 

may  write  {3.  39b)  in  the  forrh  corresponding  to  (3.  28b)  and  then  divide  by 
(3.39a).  The  result  is  ■ 


^+2-  ^  ^  V  (?> +'X‘*^+2.)  • 


Only  the  last  term  involves  gf  and  therefore  Oq*  It  is  easy  to  show  that  as 
gf  —  ow  term  is  of  order  g£  -  g^  for  g^  0  and  g£'^  for  g^  =  0. 

Tnus  the  last  term  vanishes  as  CTq  -*  0  and  we  recover  the  end  wall  result 
Eq.  ^{2.  19).'. 


For  the  second  case,  where  the  momentum  layer  is  thicker, 
Eqs.  (3.  35)  and  (3,  36)  become 


I 


(3, 4 1  a) 


m 


T'’  -I^  t’  -I 

<X.co4>l  -S-UJ-V2.  J- 


3  t5 


•C04>t 


o 


Co+J 

li-l’ 


-U>  J-Ui+t 
to-rl 


{3.  41b) 


where 


T^-  f 

Xv,  - 


(3.  42) 


For  a  given  set  of  values  g,^,  u>,  Cq  and  j(  =0  or  1),  one  must  solve 
either  Eqs.  (3.39)  for  gf  and  (d(i2/dq)w  or  Eqs.  (3.  41)  for  N  and  (d62/dq ).j(,. 
By  their  definitions,  we  must  have  gw  <  gf  <1  and  Nf  > 0,  and  this  restric¬ 
tion  has  in  all  Cases  so  far  calculated  been  sufficient  to  provide  a  unique 
solution.  The  relative  size  of  the  momentum  and  energy  layers  is  therefore 
determined  after  the  solution  is  obtained. 

One  feature  of  the  solutions  immediately  apparent  is  that  j  and  cTq 
(the  pressure  gradient  parameter  and  Prandtl  number)  appear  only  in  com¬ 
bination  and  never  separately.  This  is  a  direct  consequence  of  using  the 
simplified  condition  (3.26)  to  determine  A.  If  the  actual  momentunri  integral 
were  used,  and  j  would  also  appear  separately.  To  the  author's  knowledge, 
the  fact  that  j  and  Oq  can  be  combined  into  a  single  parameter  has  not  previ¬ 
ously  been  suggested,  and  the  only  exact  calculations  extensive  enough  to  test 
it  are  the  =  constant  ones  of  Squire,  ^  Subulkin,  ®  and  Cohen  and 

Reshotko.  From  the  results  of  these  references,  and  some  extensions 

of  them,  (di^2/'^T)w/(l  -  gw)  is  plotted  in  Fig.  2  against  (1  4  j)~  icrQ-l  for 
gw  =  P,  1.  For  each  gw,  the  correlation  is  seen  to  be  quite  good,  thus 
lending  some  support  to  the  present  approximate  theory  even  before  any 
calculations  are  rnade. 
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In  fee  incbmpressible  limit  i  »  which  can  be  deriyed  from  the 

above  equations,  since  the  value  of  does  not  enter  the  problem  in  this 
limit/;{3.  43a) -becomes  ■  - 


5  “*  ‘*‘1 


i  20*J)^  \  |  _ 


{3.45a) 


and  when  this  is  combined  with  the  limit  of  (3.  43b),  a  qua:dratic  equation  re¬ 
sults  whose  relevant  solution  is 


,  .  --i-  4-  P-L  +  -L-.  -  i  1 

:  '  I  ji_ 

(o 


(3.45bj 


Equations  {3.45)  provide  an  explicit  formula  for  the  heat  transfer  in  the 
Lncompressible,  low  Prandtl  number  case.  From  Eq.  (3.45b)  and  its 
derivative,  one  can  show  that  the  required  physical  limitation  that 
0  <  (gf  -\g..v)/(i  ■  gw)  -  i  (which  required  elimination  of  the  other  root  of 
the  quadratic)  is  satisfied  only  for  0<  (1  +  j)  Qq  <  2,  giving  the  range  of 
Prandtl  number  for  which  Eqs.  (3.45)  are  applic'ble.  For  larger  values, 
the  momentum  layer  is  thicker  than  the  energy  layer  in  the  present  theory, 
and  Eqs.  (3.44)  can  be  written 


1  _  > 


(3.46a) 


g  —  1 : 
®w 
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(3.  46  b) 


Calculation  of  the  heat  transfer  in  this  incompressible^  high  Prandtl  humber 
case  requires  solution  of  the  simple  quartic  equation  (3.  46b).^^  '  : 


By  evar.iining  the  left  side  of  (3.  46b)  and  its  deriyative,  it  can  be  verified  that 
th.e  required  limitation  N£/(  1  -  :g^v)  -  ^  4  s  olution  only  for 

(1  4^  j)crji  3.  One  can  also  see  that  at  the  critic .al  value  ,(1  +  jlap  =  2,  the 
left  side  of  {3.  46b)  changes  sigh,  thus  yerifving  in  this  simple  Case  the  niimer- 
ical  criterioh  to  be  described  below  for  deciding  which  set  of  equations  will 
provide  the  solution  for  a  given  value  of  (1  -r  j)cro* 

We  have  , already  pointed  Out  th«at  in  the  liinit  Oq -*  0,  Eqs .  (3,  39)  re¬ 
duce  to  the  proper  results,  namely  th.e  end  wall  case  considered  previously. 
The  opposite  liniit  (Jq  “  can  also  be  de rived  easily.  The  appropriate  i 
starting  point  is  EqS.  (3.  41),  since  the  niomehtum  layer  will  be  thicker  than 
the  energy  layer.  EqS.  (3.  41 )  indicate  two  possibilities  as  CTq  —  »  .  One  is 
that  N£  remain  finite,  in  which  case  (do^  '^d r)  )._y  becomes  zero  and  the  quantity 
in  curly  brackets  in  {3,  41b)  nvust  vanish..  In  the  limits  =  1  and  —  1 , 
this  leads  to  X£  <  0,  an  unacceptable  solutioh;  so  it  probably  is  unacceptable 
:  in  the  gene ral  case  als  o.  The  Othe r  pos sibility  is  th.e  more  r eas  onable  one , 
that  the  nioir.entum  laye r  thickne.ss  represented  by  X£  approaches  infinity. 

V.'ith  this  assumption,  only  thie  X£  term  on  the  right  of  (3.  41a)  and  the  first 
term  and  Xf  term  in  (3.  41b)  are  import.ant  and  a  simple  calculation  gives 


■  This  is  an  asymptotic  e.'ipressiOn  for  the  heat  transfer  rate  n-arameter  in  the 
general  case.  It  is  well  known  that  as  Oq”"  the  Prandtl  number  depend¬ 
ence  of  the  temperature  gradient  is  Equation  (3.47)  is  in  agreement 

with  this  when  we  recall  that  the  present  definition  of  rj  differs  from  the  usual 
definition  and  to  convert  to  the  latter,  one  rriust  imagine  the  right  side  of 
(3.  47)  multiplied  by  .yo^. 

The  cases  -  =  1  and  a...  —  1  both  lead  to  the  same  expression 

C^_CO  :  , 

One  can  argue  that  both  these  cases  should  have  the  same  limit.  They  haye 
the  same  energy  equation  and  differ  only  Ln  the  factor  g  instead  of  unity  in 
the  pressure  gradient  term  in  the  momentum  equation.  But  in  this  limit,  the 
momentum  layer  is  so  much  thicker  than  the  energy  layer  that  it  sees  a  con¬ 
stant  temperature  flow  at  the  external  temperature,  so  g  =1.  Thus  the 
—  T and  -  =  1  cases  satisfy  the  same  equations  and  so  have  the  same  ex¬ 
pression  for  heat  transfer.  For  more  general  cases,  where  appears  ex¬ 
plicitly  in  the  energy  equation,  one  would  expect  both  -  and  gw  to  appear  in 

the  limiting  value,  as:  they  do  in  Eq.  (3.  47).  . 

'  -Z5  - ’  ,■'■' 


Merk^^  has  given  an  exact  asymptotic  expression  for  the  case 
which  should  also  be  valid  for  the  case  ^  =  1,  according  to  the  argument  just 
given.  In  the  present  notation.  Me rk's  result  is 


^  j_  -I- 
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The  last  bracket  contains  the  wall  shear  stress  parameter,  which  is  inde¬ 
pendent  of  Oq  in  this  case.  Comparison  between  (3.48)  and  f3.  49)  will  be 
given  later  .  The  author  has  been  unable  to  find  a  comparable  exact  formula 
for  w  /  1,  g^T^  1,  so  Eq.  (3.  47)  represents  an  otherwise  unavailable  expres¬ 
sion  for  the  case  of  cooled  walls  with  pk/cp  not  constant. 

Numerical  Solutions  and  Comparisons  with  Exact  Solutions 

We  will  now  proceed  to  the  actual  numerical  solution  of  Eqs.  (3.  39) 
and  (3.  41).  This  is  best  done  on  a  digital  computer,  since  Eqs.  (3.39)  are 
transcendental  in  gf,  while  Eqs.  (3.41)  can  be  reduced  to  a  quartic  in 
(dd2/dTi)w  or  Nf.  The  advantage  in  using  the  approximate  method  over 
actually  integrating  the  exact  differential  equations  (3.  38)  is  in  the  consider¬ 
able  saving  in  computing  time.  Tbe  approximate  method  requires  only  the 
iterative  solution  of  transcendental  or  quartic  equations,  while  the  exact 
solution  requires  integration  of  a  two-point  boundary-value  problem,  with  a 
number  of  integrations  to  q  »  1  usually  necessary  before  the  outer  boundary 
condition  is  satisfied. 

The  procedure  used  to  obtain  a  solution  was  the  following.  The  value 
gf  -  I  was  used  in  (3.  39a)  to  obtain  a  value  of  {dd2/‘^q)w>  was  inserted 
into  the  left  side  of  (3.  39b)  together  with  g£  =  1.  If  the  left  side  of  (3.  39b) 
was  positive  with  these  values,  a  new  guess  for  g£  <  1  was  made,  and  the 
procedure  repeated  until  (3.  39b)  was  satisfied,  which  always  occurred  for 
gvi^  <  g£  <  1.  If  the  left  side  of  (3,  39b)  was  negative  for  g£  >1,  attention  was 
shifted  to  Eqs.  (3.41),  which  were  then  solved  by  iteration,  starting  at 
Nf  =  0  and  searching  through  increasing  values  of  N£.  This  method  has  pro¬ 
duced  a  solution  in  all  cases  attempted  so  far  and  was  verified  analytically 
in  the  discussion  of  the  limiting  case  g^  —  1  given  above. 

To  compare  the  solutions  of  the  present  approximate  theory  with  ex¬ 
act  solutions  of  Eqs,  (3.  38),  use  may  again  be  made  of  the  results  of  Refs. 

8,  9i  and  10  as  well  as  recent  work  of  Bade.  Such  a  conniparison  is  shown 
in  Table  III.  Where  no  reference  is  given,  the  integration  of  the  exact  equa¬ 
tions  was  done  by  the  author. 
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Table  III  shows  the  same  striking  result  as  Tables  I  and  II,  namely 
the  small  percentage  variation  in  the  ratio  of  exact  to  approximate  values  of 
(d(J2/‘i’l)w  restricts  oneself  to  values  of  (1  +  j)CTo  >  0*  i»  the  varia¬ 

tion  is  from  1.  02  to  1.  10,  or  a  +  4%  variation  about  1.  06.  Thus  a  suitable 
correction  factor  in  this  range  Is  1.06.  To  make  this  fair  smoothly  into  the 
factor  i.  13  suggested  for  the  end  wall  case  (cTo  — ■  0),  one  can  add  a  term 
which  is  negligible  for  (1  +  j)CTrt  >  0.  1  but  approaches  0.  07  as  ^  0.  A 
suitable  correction  factor  which  fullfills  these  req^uiremehts  is  1. 06  +  0.  07 
exp  [-10(1  +  j)cro]  •  Use  of  this  factor  changes  the  heat  transfer  relation 
■|3.l7):tO;: 


(3.50) 


which  has  a  maximum  error  of  +  4%  for  all  the  cases  in  Table  III.  Use  of  this 
formula  with  the  present  approxTmate  method  of  solution  of  the  stagnation- 
point  boundary-layer  equations  yields  an  accurate  prediction  of  heat  transfer 
rate  for  power  law  fluid  properties  and  constant  Prandtl  number.  The  neces¬ 
sary  computing  time  is  several  factors  of  ten  less  than  that  required  to  solve 
the  exact  two-point  boundary- value  problem. 
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SECTION  4 
CONCLUSIONS 


A  simple  integral  method  has  been  suggested  for  calculation  of  heat 
transfer  rate  from  similarity  boundary  layer  equations.  The  method  has  been 
applied  to  the  end  wall  of  a  shock  tube  behind  a  reflected  shock  (  Rayleigh-type 
problem)  for  a  perfect  gas  and  to  the  stagnation  point  for  a  gas  in  thermody¬ 
namic  equilibrium  with  Lewis  number  unity.  In  both  cases  explicit  expres¬ 
sions  were  derived  when  the  thermal  conductivity,  density,  and  specific  heat 
were  expressed  as  powers  of  the  temperature  or  enthalpy,  and  the  Prandtl 
number  was  constant. 

For  the  end  wall  geometry,  a  simple  analytical  expression  for  the 
heat  transfer  was  obtained  and  compared  vvith  exact  calculations.  A  constant 
correction  factor  brought  the  formula  into  agreement  with  the  exact  results 
within  +  3%  for  a  large  range  of  powers  and  wall  to  sLream  temperature 
ratios.  The  expression  for  the  heat  transfer  rate  is 


where  (dp j/dq is  obtained  from  Eq.  (2.  l6a).  In  the  case  T“,  Cp  con¬ 
stant,  Eq.  (2.  19)  gives  an  explicit  formula  for  (d<^)j/dq)^  already  obtained 
by  Jepson.  ^  In  the  case  where  k  is  proportional  to  two  different  powers 
of  T  in  different  temperature  rarsges,  (d(;ii/dq)^y  is  found  from  Eq»  (2,  23). 

The  same  correction  factor  (1.  13)  holds  for  both  cases  and  presumably  for 
any  linear  combination  of  power  laws. 

For  the  stagnation  point  geometry,  two  pair  of  transcendental  equa¬ 
tions  (3.39)  and  (3,41)  were  derived,  one  pair  of  which  in  any  given  case 
provides  the  solution  for  the  heat  transfer  parameter  and  relative  size  of 
the  momentum  and  energy  layers.  The  two  pair  of  equations  must  be  con¬ 
sidered  together,  since  one  is  appropriate  when  the  momentum  layer  is 
thinner  than  the  energy  layer,  the  other  appropriate  in  the  reverse  situation. 
Machine  solution  of  these  equations  by  trial  and  error  is  many  times  faster 
than  solution  of  the  exact  differential  equations.  A  feature  of  the  approximate 
equations  is  that  the  Prandtl  number  and  pressure  gradient  parameter  only 
appear  together  as  one  parameter.  In  the  only  case  known  to  the  author  where 
exact  calculations  cover  an  extensive  enough  raiige  of  these  two  parameters 
to  permit  testing  this  result,  namely  the  case  pk  /cp  =  constant,  the  single 


-29- 


parameter  correlated  the  heat  transfer  rate  within  a  few  percent. 

When  the  results  of  the  approximate  stagnation  point  calculations  for 
power  law  fluid  properties  and  constant  Prandtl  number  were  compared  to 
exact  values,  a  correction  factor  was  found  which  brought  agreement  within 
t  4%.  The  heat  transfer  rate  given  by  the  aooroximate  theory  is 


>2:^^  \d%4 


Jh 


fr  >  -fo(i+iK 

\  |.o^>+o.ole 


where  (d(i»2/dr])^,,  is  obtained  by  solving  the  appropriate  pair  of  Eqs.  (3.  39) 
or  (3.41). 
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